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We studied the finite temperature phase transition of SU(2) gauge theory on 
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our fiux data with predictions of Michael sum rules, we studied the behavior of string 
tension with temperature. Our data agree well with the generalized sum rules. 
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I. INTRODUCTION 



The finite temperature phase transition of QCD [1] has been studied extensively in 
lattice gauge theory (LGT). It has been shown that at low temperatures QCD is in the 
confined phase, quarks are confined by a linear interquark potential, V{r) ~ Kr, with k 
the string tension. However, at sufficiently high temperatures the system is in the uncon- 
fined phase. Quarks are unconfined because the interquark potential is a screened Coulomb 
potential [2-4], V{r) - e'^Vr. 

It is believed that quark confinement is due to the string formation in the confined phase, 
that is, the color field between quarks forms a string-hke fiux tube. In LGT one can probe 
the spacial distribution of the energy density between a qq pair by using a local operator 
(e.g., plaquette). Some strong numerical evidence for string formation has been found [5-8]. 
Here we measured the fiux data at various temperatures. The color fiux energy and action of 
the qq pair are related to the interquark potential energy V{r) by the Michael sum rules [9]. 
In this paper our aim is to check these sum rules with our finite temperature fiux data. 

The Michael sum rule for energy is the LGT version of the relation in classical electro- 
statics in which the work to assemble a charge distribution is equal to the integral over all 
space of the electrostatic energy density. The Michael sum rule for energy and a similar 
relation for action were tested on large lattices by Haymaker and Wosiek [10]. Except for 
an unresolved discrepency arising in the determination of self energies, both sum rules were 
satisfied by the fiux data. In this paper we rederive the sum rules for the case of finite 
temperature and check them with our fiux data. 

In the following we shall only consider pure SU(2) gauge theory. To estimate the pre- 
diction of the Michael sum rules at finite temperatures we need to know the lattice scaling 
relation a(/?) and the behavior of the string tension with temperature, k{T). For this pur- 
pose we adopted the lattice asymptotic scaling relation a(/9) [11] to the non-perturbative 
region by allowing the scale constant to vary slowly with (3. Also the string tension was 
measured at finite temperature, and it was found that the string tension data agree very 
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well with the fitting function, 

K{T)^Ko{l-^r (for T<T,). (1) 

^ c 

This paper is organized as follows. In Sec. II we review some basic properties of the finite 
temperature phase transition for SU(2) gauge theory. The scaling relation a{P) is studied, 
and the flux measurement techniques in LGT is presented. In Sec. Ill we describe our 
Monte Carlo calculations, and present some numerical results for global quantities, then we 
analyze the behavior of string tension k with temperature. In Sec. IV we present a complete 
derivation for Michael sum rules and generalize them to include finite size effects, then we 
give a detailed comparison of these sum rules with the flux data. Finally Sec. V presents 
the summary and conclusions. 

II. THEORETICAL BACKGROUND 

In this section we present the fundamental ideas of LGT in studying the finite temper- 
ature phase transition of QCD, and emphasize aspects that are relevant for our numerical 
investigations. 

A. SU(2) Finite Temperature Phase Transition 

In our study of SU(2) LGT we used the standard Wilson action [11] 

5(t/)=/?^(l-lTr^p), (2) 
p ^ 

where /? = A^j with g to be the lattice coupling constant, and C/p is the product of link 
variables around a plaquette, 

Up{n) = U,{n)U,{n + iAU;\n + v)U-\n), (3) 

where the indices /x, u represent the orientation of the plaquette, n the position of the 
plaquette on the lattice. 
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On a four-dimensional Euclidean lattice of the size, Nt x x Nz, the temperature T is 
related to the temporal extent of the lattice [12]. If one chooses Nt as the time direction, 
then one can define the temperature as. 



where a is the lattice spacing, which is a function of (3. 

To study the finite temperature phase transition a convenient order parameter is provided 
by the Polyakov loop closed in the time direction [13], 



where C/t(r, r) is the link variable oriented in A^^ direction at the site (r, r). It is well-known 
for pure SU(A'") gauge theory that the expectation value of the Polyakov loop (P) plays the 
role of an order parameter [13] in the finite temperature phase transition, that is, in the 
infinite volume limit one has 

in the confined phase, T < Tc, {P) = 0; 

in the unconfined phase, T > (P) 7^ 0; 
because of the spontaneous breaking of the global Zn symmetry [14] . Here Tc is the transition 
temperature. 

The expectation value (P) is associated with the free energy Fq of an isolated quark q 
in the infinite volume limit [14], 



where Lt — Nta, is the inverse temperature. Further, the correlation function (P(O)P^(r)) 
is related to the potential energy V{r) of a qq pair in the infinite volume limit. 



1 JVt 

P(r) = -Trn^*(r,r), 



(5) 



T=l 




(6) 



(P(O)P^(r)) ~ canst. ex.p[-LtV{r)]. 



(7) 



In the confined phase, T < Eq. (7) becomes, for large r [2], 



lim (P(O)P^(r)) ~ exp[-LtKr], 
= {Pf = 0, 



(8) 
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because of (P) = in the confined phase, where k is the string tension. Eq. (8) imphes that 
the qq potential in the confined phase is hnear, e.g., V{r) ~ kt. On the other hand, in the 
unconfined phase, T > Tc, and in the infinite volume limit it is found that for SU(2) gauge 
theory [2], 

3 

lim (P(O)pt(r)) = const.[l + L^— — exp(-mr)] 

= {P? + 0, (9) 

for (P) 7^ in the unconfined phase. Eq. (9) imphes that the qq potential in the unconfined 
phase is a screened Coulomb potential, e.g., V{r) ~ with the Debye screening 
length. To the lowest order of perturbation theory one has [2] 

= y^T^ (10) 

Ref. [3] gives a detailed discussion about the qq potential in the unconfined phase. We shall 
not discuss it any more in the following. 

The transition temperature can be determined by Monte Carlo calculations of LGT. 
For SU(2) one recent result is [15,16], 

/3c = 2.2985 ± 0.0006. (for Nt = 4) (11) 

This transition point /3c was calculated on large lattices (e.g., 4 x 26^), on which finite size 
effects are small, the results can be considered to be the transition point in the infinite 
volume limit. 

To obtain the transition temperature Tc in physical units, one must know the scaling 
relation between the lattice spacing a and /3. If we assume the scaling relation a(/3) is given 
by the asymptotic relation for SU(2) [11], 

= AZ'(^vr2/5)^ea;p[-A^V] = A-i/(/3), (12) 

where is the lattice scaling constant, then one can calculate the transition temperature 
Tc in physical units from Eqs. (4), (11) and (12), 
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Ntfm' 



= 42.11 ±0.06. 



(13) 



This is consistent with the results of Refs. [17,22]. From Eqs. (12) and (13) one can also 
see that P < Pc implies T < Tc, and P > Pc corresponds to T > Tc for fixed Nf. Now we 
proceed to study the lattice scaling constant A^. 



We follow Refs. [19,20] in choosing a scale such that in the limit of zero temperature and 
infinite volumes the string tension is 



which was determined from the Regge trajectory of experimental data. Here for definiteness 
we use the data from the real world to discuss SU(2) gauge theory. 

Prom Eq. (14) and the string tension data of Monte Carlo calculations in Refs. [8,19] 
one can extract a relation between the lattice spacing a and lattice coupling constant p. 
The result is listed in Table I, which is consistent with the similar result of Ref. [10]. In 
Table I the values of were calculated from the asymptotic scahng relation in Eq. (12) 
by substituting the corresponding values of a and P in this table. One can see that the 
values of A^^ changes slowly with P in the region, 2.22 < P < 2.5. This implies that the 
perturbative asymptotic scaling relation in Eq. (12) is not exactly valid in this region. To 
study the non-perturbative physics and obtain the relation a{P) for other P values in this 
region, we use the function in Eq. (12) to fit the data of Table I, with A^^ considered as a 
function of p. For simplicity, we chose a quadratic function to fit the A£^ values. The result 
is 



B. The Scaling Relation a{P) 



^ = 0.44 GeV, 



(14) 



a{P)=Al\P)f{P), 
with Al\p)^di + d2p + dsp^ (fm), 
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(15) 



where the coefficients di = 59.37 ± 0.86, = -5.96 ± 0.39 and = -3.28 ± 0.17, and the 
function f{P) is given by Eq. (12). 

Using the scahng relation of Eq. (15) we can transform the Monte Carlo data from lattice 
units to physical units. Now the transition temperature Tc in Eq. (13) can be expressed 
explicitly in physical units, i.e., 

Tc = 1-487 ± 0.140 (1/fm). (16) 

Also, by using the scaling relation one can estimate the transition point Pc on lattices of any 
temporal size Nt, if we assume the transition temperature Tc in Eq. (16) is independent of 
the lattice size Nt because it is a physical observable. For example, the estimated values of 
Pc for Nt — 6 and 8 are 

Pc ^ 2.42 (for TV* = 6), 

pc ^ 2.50 (for Nt = 8). (17) 

In Sec. IV we shall use the scaling relation of Eq. (15) to discuss the predictions of Michael 
sum rules. In the following we proceed to describe the methods used in measuring the qq 
flux distributions. 



C. qq Flux Distributions 



To measure the flux distributions of a qq pair, we calculated the quantity [6,8], 



/^,(r,x) = — 

til 



(P(0)Pt(r)n^,) 



- M 



(18) 



(P(O)Pt(r)) 

where the Polyakov loop P(r) is defined in Eq. (5), and D^,^ = |Tr(C/p) is the plaquette 
variable with Up defined in Eq. (3). The plaquette has 6 different orientations, (/x, i') = {2, 3), 
(1,3), (1,2), (1,4), (2,4), (3,4). 

To reduce the fiuctuations of the quantity P(0)P'^(r)n, in practical calculations we mea- 
sure the quantity [6] 
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(P(0)Pt(r)n^,(x)) - (P(0)Pt(r)n^,(x^)) 



(19) 



(P(O)Pt(r)) 

as the flux distribution instead of Eq. (18), wliere tfie reference point x^j was chosen to be 
far from the qq sources. This replacement does not change the measured average due to the 
cluster decomposition theorem. The six components of /^^ in Eq. (19) correspond to the 
components of the chromo-electric and chromo- magnetic fields {S,B) in Minkowski space, 
i.e., 

- li-13l -Bl -Bl £l El El). (20) 
The total electric and magnetic energy densities are defined as 

Pr..^\{B') = \m) + {Bl) + {Bl)]. (21) 

To improve the statistical accuracy we used the multihit techniques [6] , that is, we make 
the replacement for Polyakov loops in Eq. (19), 

1 Nt _ 

P(r)^P(r) = -Trn^*(r,r), (22) 

^ T=l 

where t/t(r, r) is given by 

- _ !UteM(^l^{UtXl + h.c.)]dUt 
* !eKY>[l3T^{UtXl + h.c.)]dUt ' 

with X} to be the neighborhood of C/j, i.e., UtX} = J] I^Pj the sum extends over all plaquettes 



containing Ut- And A — ^det{Xt), h and I2 are the modified Bessel functions. This 
technique is limited to the situations where the neighboring links Xf of Polyakov loops P 
do not overlap with the plaquette O^^, in Eq. (19). 



8 



III. NUMERICAL RESULTS FOR GLOBAL QUANTITIES 



In this section we present the numerical results and analyses for global quantities in 
our Monte Carlo calculations. We follow Ref. [17] to choose our lattices of the geometry 
Nt X Ng X Nz, with Nt <^ Ng <^ N^. If Nt is chosen as the time direction, we can simulate 
SU(2) gauge theory at finite temperatures {T — l/Nto) in large volumes {V — NgN^a^). 
Typically, we choose A^t = 4 and 6, Ng = 7, 9, 11, and Nz = 65 in most cases and = 37 
in a few cases with A^^ — 6. The lattice coupling constant /? is in the range 2.25 < {3 < 2.40 
for Nf — 4, and 2.30 < (3 < 2.50 for Nt — 6. We updated the lattice configurations 
by using the standard Metropohs algorithm alternated with overrelaxation method. We 
typically thermalized for 4000 sweeps, and made one measurement every 10 sweeps. The 
total number of measurements for each data is about 2000. The actual number may vary by 
a small amount in each case. The calculations were done on LSU's IBM 3090 mainframe. 

A. Measurements of The Order Parameter 

and {P{{))P^{z)) 

As we discussed in Sec. II A the expectation (P) plays the role of an order parameter 
in the infinite volume limit. However, on a finite lattice this quantity is always zero if 
the computation time is taken to be infinite, because the system would flip between two 
ordered states, and the values of the Polyakov loop P would flip sign after some iterations. 
Therefore, we take the expectation value of the modulus of the Polyakov loop (|P|) as the 
"order parameter" on finite lattices [13,14]. 

The measured data of (|P|) are plotted in Fig. 1, which shows a rapid increase of (|P|) 
at /3 a; 2.30 for Nt — 4, and another at /? fti 2.42 for Nt = 6. This implies that a phase 
transition occurs at Pc ~ 2.30 for Nt — 4, or Pc ~ 2.42 for Nt — 6, in the infinite volume 
hmit, which are consistent with Eqs. (11) and (17). For cases of A^t = 4 the data measured 
on the lattice 4 x 9^ x 65 agree very well with those from 4 x 11^ x 65, both approach the 
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infinite volume limit. However, for cases oi Nt — 6 the data from the lattice 6 x 7^ x 65 
have some discrepencies with the data from 6 x 11^ x 37. This may be due to finite- volume 
effects, as discussed in Ref. [15], because in these cases Ns (or N^) are not large enough to 
approach the infinite volume limit. 

We also measured the correlations of two Polyakov loops, (P(O)P^(z)), with the sepa- 
ration z along Nz direction. Some typical results are shown in Figs. 2 and 3, both were 
drawn as a logarithmic plot. Fig. 2 shows (P(O)P^(r)) versus r in the confined phase, from 
this figure one can see that the correlation (P(O)P^(r)) — > as r — > oo, as predicted by 
Eq. (8). Also, the linear behavior of the quantity ln(P(0)P^(r)) with r indicates that the 
qq potential V{r) is linear in the confined phase. And the slope of the straight lines corre- 
sponds to the string tension k. This figure shows that the string tension decreases with /3 
(or temperature). Fig. 3 is similar to Fig. 2, but is in the unconfined phase. It shows that 
the correlation (P(O)P^(r)) 7^ as r — > oo, as predicted by Eq. (9). 

B. Measurements of String Tension k 

The string tension k can be extracted from the correlation (P(O)P^(r)) according to 
Eq. (8). In our measurements we choose the correlation along the longest extent N^, i.e., 
r = re^, and follow Ref. [17] to extract the string tension data from the relation, 

{P{0)P\z))c = const.{e-^*^' + e-^*«(-^^-^)), (24) 

where on the L.H.S. of this equation the connected correlation {P(0)P'^(2;))c is defined as, 
{P{0)P\z)) — (P)^. So the dependence on the coordinates t, x, y can be suppressed, and one 
can extract the effective string tension. On the R.H.S. the new term q-^Ml^-z) accounts 
for the finite size effects along the correlation direction A^^, with — N^a. 

The raw string tension data are given in Table 11. Some data were measured in the 
transition region, e.g., /3 = 2.30 for (Nf = 4). In this case we assume Eq. (24) be still vahd. 
Since our lattices have the geometry, ^ Ng ^ Nt, we find that finite size effects are 
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small as Ng/Nt > 2 in the limit of — > oo. The detailed analysis about the finite size 
effects is presented elsewhere [5]. This conclusion also agree with the results of Refs. [17,18]. 
Therefore, we expect that the data measured on lattices 4 x 9^ x 65 and 4 x 11^ x 65 can be 
viewed as the string tension in the large volume limit. However, the data from the lattice 
6 x 7^ x 65 might have some large finite size effects because A^^ = 7 is not large enough for 
Nt = 6. For data measured on the lattice 6 x 11^ x 37 the transverse size A^^ = 11 is large, 
but some large finite size effects may be caused by the relatively small longitudinal size Nz- 



C. Behaviour of k With Temperature 

In Sec. IV when we discuss Michael sum rules and analyze the finite temperature effects 
on the fiux data, we need to know the relation between string tension k and temperature. 
In this part we proceed to analyze the string tension data. 

Using the scaling relation in Eq. (15) we can transform the string tension data in Table II 
from lattice units to physical units. The results are shown in Table III. Prom this table one 
can see that the string tension k decreases with temperature T (or /3) in the confined phase. 
If we assume that «; is a continous function of T (or one has 

^<0 (forT<T,), 
or — <0 (for/3</3,). (25) 

This behavior is also confirmed in other recent papers [21,22]. 

Since the string tension k decreases with temperature T in the confined phase and is 
expected to vanish at the transition point Tc, for simplicity, we then use the simple function 
with power behavior in Eq. (1) to fit the string tension near Tc, that is, 

k{T)^Ko{1-^Y {iorT<n). (26) 

c 

This is a simphfied form of the fitting function used in Ref. [22]. Here the constant kq 
represents the string tension k in the limit of zero temperature and infinite volumes, which 
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can be given by Eq. (14), i.e., kq ~ 0.981 GeV/fm. The transition temperature Tc is given 
by Eq. (16) in physical units. 

As we discussed above, the string tension data measured on lattices 6 x 7^ x 65 and 
6 X 11^ X 37 may have large finite size effects. By viewing the data in Table III we find that 
the string tension data from the lattice 6 x 7^ x 65 are consistent with the data from lattices 
4 X 9^ X 65 and 4 x 11^ x 65, on which finite size effects can be neglected, as we discussed 
before. However, data from the lattice 6 x 11^ x 37 have large discrepancies with other data. 
This suggests that finite size effects on string tension are sensitive to the longitudinal size 
Nz, but not as much to the transverse size Ns- 

We then use the string tension data from lattices of small finite size effects (i.e., 4 x 11^ x 65 
and 6 X 7^ X 65), to fit the function in Eq. (26). The only fitting parameter is the exponential 
index a, the best fit result is, 

a = 0.35 ± 0.04. (27) 

Fig. 4 shows this fitting function and the measured string tension data of small finite size 
effects. In this figure one can see that a few data point very close to the transition point 
Tc have large discrepancies from the fitting function, because fiuctuations are large in the 
transition region. We then excluded these data at, T pa 1.48 (1/fm), from our fitting 
process. Also all other data above Tc were excluded because they would correspond to 
complex values for the fitting function. Below the transition region our data agree very well 
with the conjectured behavior in Eq. (26). 

To get the relation between k and (3, we also fit the string tension data with the following 
function similar to Eq. (26), 

«(/3) = «o(l - f )' (for/3</3c). (28) 

The transition point Pc = 2.2985 ± 0.0006 for A^* = 4 can be chosen from Eq. (11), and 
/3c = 2.42 ± 0.01 for Nt = 6 from Eq. (17). We find that the best fit results for the 
exponential index S are 
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5 = 0.22 ± 0.03 



(for TV, = 4), 



6 = 0.14 ±0.04 



(for Nt = 6). 



(29) 



In Fig. 5 we plot the result for the case of Nf = 4. Also, in this figure the data above 
the transition point, /3c — 2.2985, were excluded from our fitting process. Again our data 
agree well with the fitting function. In Sec. IV B we shall use the relation k(/3) to check the 
Michael sum rules. 



To study the detailed properties of the finite temperature phase transition, we measured 
the qq flux distributions at various temperatures using the technique discussed in Sec. subsec: 
mfd. The flux energies (action) arc related to the potential energy V{r) of the qq pair by 
Michael sum rules, which were first derived by C. Michael [9]. In the derivation of these sum 
rules some scaling relations on asymmetric lattices are used, which were studied extensively 
by A. Hasenfratz and P. Hasenfratz [24], and F. Karsch [25]. The original work of C. Michael 
is somewhat terse and drawing on other references. In order to generalize it we feel it is 
helpful to present a thorough derivation, including the other references. In this section we 
shall give a complete derivation of Michael sum rules, and present our supplements and 
generalizations. One result is delegated to the Appendix. Then we compare our flux data 
with the predictions of these sum rules. 



For a static qq system at the spatial separation r the Michael sum rules states that, in 
the limit of zero temperature, infinite volume and infinitesimal lattice spacing (a — > 0), the 
electric and magnetic flux energies of the system in Minkowski space are [9], 



IV. MICHAEL SUM RULES AND FLUX DISTRIBUTIONS 



A. Review and Generalizations of Michael Sum Rules 



(30) 
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where the sums extends over the whole space around the qq pair, V{r) is the potential 
energy and //a is the self energy of the qq sources. Eei{r) and £'ma(^) denote the electric 
and magnetic parts of the flux energy, and (B^) are given by Eq. (21). 



1. The Action Sum Rule 

To derive Eqs. (30) and (31), it is convenient to work in the transfer matrix formalism 
of LGT. We choose the temporal gauge [23], that is, all temporal links are trivial, C/t(n) = 
gigaA^(n) _ ^ Q^g remark about the temporal gauge of a finite lattice is that along a time 
axis one can not choose all temporal links trivial, Uf = I, because of the restriction that 
a Polyakov loop in the time direction is gauge invariant, which has values other than the 
trivial one, P{f) — 1. Therefore, on a finite lattice one can choose all temporal links to be 
trivial, except one link on each time axis, which can be chosen far from the operators under 
considerations. 

Let's consider a Wilson loop W of the temporal size na and the spatial size r, as shown 
in Fig. 6. The time directed pathos represent the qq sources, the space directed pathes Pr{0) 
and Prin) create and annihilate the static qq pair from and to the vacuum. In the temporal 
gauge the expectation of Wilson loop becomes, 

(W) = {Pr{0)Pr{n)), (32) 

where the expectation is evaluated in the partition function form. 



Jd[U]e-'^'^'W 

7W¥- 



with pS' — S, the action given by Eq. (2), i.e., S' — Z](l — □), where □ is defined in Eq. (18), 
the sum is over all plaquettes on the lattice. 

Using the transfer matrix approach [23,7] and assuming a discrete spectrum for the lowest 
eigenstates of the transfer matrix, one can evaluate Eq. (32) as 
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1 , 

1 



{W) = ^Tr(n(0)r^^-P,(n)r^-") 



= ^tA^E^.^^^.^"W<*""' (34) 

where Tqq is the transfer matrix projected into the qq sector of the Hilbert space [7], T is the 
transfer matrix for the vacuum. r) and |//) are the eigenstates of Tqq and T respectively, 
and Ay(r), A^^ are the corresponding eigenvalues, with Ajy(r) = e""^''^'') and A^ = e~°'^^. The 
coefficient d^i, = {jj,\Pr\i',r). In the limit of A^j — > oo, the partition function Z = Xq \ and 
in Eq. (34) the dominant contributions correspond to /i — 0, that is, 

W = E^Lf^)" (asTV.^oo). (35) 



Ao 

If the temporal size n of the Wilson loop W is very large {n oo), the dominant term in 
Eq. (35) is given by u — 0, 

{W) = rf^o(^)" = dl,e--^(^o(r)-Eo)^ (asN^n^oc) (36) 

where the energy of the vacuum Eq is usually chosen to be zero, in the following we will 
take this choice which implies Aq = 1. 

Now we consider the /3-derivative of Eq. (36), that is, 

^ = -{WS') + {W){S') = ^^6-'^^'^^% (37) 

where wc have taken Eq = 0, and S' is defined after Eq. (33). 

Consider a plaquette □(m) outside the Wilson loop W in the time direction (i.e., < 
n < m), e.g., the plaquette Pi in Fig. 6, where we draw a plaquette with a time extension. In 
the limit of infinite large n and infinitesimal lattice spacing, a — > 0, one can neglect the time 
extension of Pi. The contribution of Pi to the L.H.S. of Eq. (37) is {Wn{m)) - {W){a{m)). 
As □(m) is far from W ( m—n ^ 1 ), one has (VFn(m)) fa (VF)(n(m)), then the contribution 
vanishes, {WD(^m)) — {W){0(^m)) — > 0. However, when □(m) is close to ( m ~ n ), one 
expects that this contribution does not vanish. Now let us show this explicitly, 
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{Wn{m)) - {W){0{m)) (forO<n<m) 

= iTr(P,(0)T^^-P,(n)T™-^^n(m)T^-™) - {W){a{m)) 

E ^0. (^) ^d,. (^) □(m) |0) 

-i:4(^)'^(0|aM|0), (38) 

where in the last step we have used Eq. (35) 

Prom Eq. (38) one can see that, as m — n » 1, the dominant term of {WD(m)) is given 
by u = 0, which would be cancelled by the product {W){0{^m)). So one has {WO[m)) — 
{W){D(m)) 0, as m - n > 1. 

However, for m fa n the term of = of the quantity, {WD(m)), in Eq. (38) is cancelled 
by the product (W) {D(^m)) . The major contribution comes from the term oi v = 1. In the 
limit of n — > oo Eq. (38) becomes, 

{WD{m)) - {W){D(m)) ( for < n< m ) 

^-uaEoir) rf^^rf^^ ^r"" (1 1 □ (m) 1 0) . (39) 



n— ►oo 



d,,[fj (l|n(m)|0) 



where in the last step we have used the fact that Aq = 1, and the dominant term is // = as 
n ^ oo. Eq. (39) implies that when the plaquette □(m) is close to the Wilson loop in the 
time direction, it gives a contribution of the order e~""^o(''), because the coefficients doo, doi 
and A7*~" are of order of unity in this case. 

Finally for the plaquette □(m) inside the Wilson loop W in the time direction (i.e., 
< m < n), e.g., the plaquette P2 in Fig. 6, its contribution to the L.H.S. of Eq. (37) is. 



{Wa{m)) - {W){D{m)) (0<m<n 
1 
Z 



= -iv(p,(o)7;-n(m)r^r"^p.(n)r^*-") - W(n(m)) 



i Vt— >o o 



V- , , /A„ r \™/A^ r A"-™, , , , 
E^OMC?.o(^j {i^,r\n{m)\u,r) 

^{W){a{m)). (40) 
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In the limit of n — > oo, and for cases that both m and n — m are large, the dominant term 
of Eq. (40) is given hy ji — v — that is, 

{WU{m)) - {W){n{m)) (0<m<n) 
-Ao(r) 



4(^)"[(0,r|n(m)|0,r) - (0|n(m)|0)] 

dU^T i^)r-o, (41) 



Ao 

where we denote (□)r-o = [(0, r|n(m)|0, r) — (0|n(m)|0)], and Eq. (36) is used. 

From the above discussion we can conclude that for plaquettes outside the Wilson loop 
W in the time direction, their contributions to the L.H.S. of Eq. (37) can be neglected as 
they are far enough from W. For the Wilson loop of large temporal size, i.e., n — >^ oo, 
the major contribution to the L.H.S. of Eq. (37) comes from plaquettes inside W. The 
contribution of one plaquette inside W is given by Eq. (41). Summing over contributions 
from all such plaquettes gives the dominant term of the L.H.S. in Eq. (37), that is, 

d{W) 



= -{WS') + {W){S') 

^4f^)"(E°)-o, (42) 



dp 

Ao 

where the sum is over all plaquettes in the spatial volume at one fixed time value m ( 
< m < n ) . The factor of n comes from summing equal contributions from plaquettes on 
each time slice. The dominant term of the R.H.S. in Eq. (37) is, as n — > oo, 

|[d^oe-"""°(^^] = -nrf^oe-'^"^"^^)^!^^. (43) 



Collecting Eqs. (37), (42) and (43) yields, 

(E °)-o, (44) 



daE^ ir) 



d/3 

where the R.H.S. can be measured by calculating the quantity, {WO)/{W) in LGT. In the 
continuum limit (a — > 0) the R.H.S. of Eq. (44) becomes, 

(E °)-o = ^ E IA{^\-o - {B\-o] = ^A, (45) 
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where £ and B are the color electric and color magnetic fields in Minkowski space, A is the 
integration of the action density over the spatial volume. 

On the L.H.S. of Eq. (44) if we write the qq colour field energy Eo{r) in the form [9], 

Eoir) = V{r) + f{i3)/a, (46) 

where V{r) is the potential energy, f/a the self-energy. Then one can obtain the following 
relation from Eq. (44) in the continuum limit, 

f3 daEo{r) _ p daV{r) (3df 

a dp a dp a dp' ^ ^ 

This is just the Michael's action sum rule, which can be obtained by subtracting Eq. (31) 
from Eq. (30) because oi A — E^i — -Bma, if we assume that V{r) scales and hence is 
independent of p. 



2. Energy Sum Rules and Generalizations 

To derive the energy sum rules in Eqs. (30) and (31) we need to study the color electric 
and color magnetic fields separately. Let us consider an asymmetric lattice with the time- 
spacing ati the spatial-spacing {— a) and the asymmetry ^ — ag/at [9]. The action for 
SU(A^) LGT becomes, 

SA^-PtY.^t-PsY.^s + const. , (48) 

where is a plaquette with a time extent and is a space-like plaquette. In the continuum 
limit (a^, at — > 0) the action must become the classical action, 

SA^\j d'x{F;,)\ (49) 

In this limit the time-like plaquette becomes 

PtU,^ PtU^^ = P.^tr^e^'ta^stF,,^ 

^ Pt{l - ^alalgl{F^,)% (50) 
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where N is for SU(A^). j — 1, 2, 3 and c is the colour index. To get the correct continuum 
action, Eq. (49), we require that in Eq. (50), 

2 2 2 3 

—ata^Qt = a^at, 
from which we obtain the relation between /St and gt- 

A ^ ^ (51) 

9t at gt 

Similarly one can obtain the relation between (5s and gs, 

p, = ?Ii-^.^4^-K (52) 

g; as gi 

In the weak coupling limit, Z^^, (5t can be expanded in terms of the coupling (5 of the 
corresponding symmetric lattice [9,24], 

2N 

M- — - P{a) + 2Ncs{0 + 0{(3-'), (53a) 

g; 

2N 

Ptr' = — = /3(a) + 2Nct{0 + 0{(3-'); (53b) 
9t 

where the coefficients Cs(^) and Ct(^) satisfy the conditions [9,25], 

c,(Ok=i = Q(Ok=i = o, 

^[csiO + ctm=i = {c's + c;)5=i = bo = 11N/A8n'. (54) 

where c'^ (cQ denotes the ^-derivative of Cg (ct)- 

With these definitions we can now study the color electric and magnetic fields of a qq 
sources. Again we consider a Wilson loop W of the size nat in the time direction and r in 
the space direction on the asymmetric lattice. We can obtain the following result in this 
case, if we repeat the steps that lead to Eq. (36), 

{W) = (i^oe-""'^''^''), (as Nt,n^oo) (55) 

where Eo{r) is the ground-state energy of the qq pair, and (W) is defined by the partition 
function formalism in Eq. (33), with the action (3S' replaced by Sa in Eq. (48). 
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By taking the derivatives of Eq. (55) with respect to (5s and (5t respectively, one can get 
the following relations similar to Eq. (44), 

datEQ{r) 



d(3t 
datEo{r) 



= (^nt),_o, (56a) 

s 

= (^n,),_o. (56b) 



where the sum J2s is over the whole spacial volume as before. The R.H.S. of Eqs. (56a) 
and (56b) become the total color electric and color magnetic energies of the qq pair respec- 
tively in the continuum limit (a^, — > 0). On the symmetric lattice (^ — 1) one has, 

lmi(i: DO.-o = ^ E la'{S\^o, (57a) 
lim(E □.).-o = E la\B')r-o; (57b) 

where S and B are the color electric and magnetic fields in Minkowski space. 

To evaluate the L.H.S. of Eqs. (56a) and (56b) we need to resort following relations 
which relate the quantities on the asymmetric lattice with the quantities on the equivalent 
symmetric lattice, 

(dmM] _i(dF l^F^ 

f dF{Ps,Pt) \ If OF ldF\ 

[-^p^),.r2[dp-pW,^^^ ^'''^ 

where F{Ps, Pt) is a function of Ps, Pt- The proof is given in the Appendix. 

Since we can write the energy Eq^v) of Eqs. (58a) and (58b) in the form, £'o(r) — 
+ f{P)/a, from Eq. (46), then the L.H.S. of Eq. (56a) becomes, 

f datEo{r) \ _ ( datV{r) \ , ( d{f/0 \ 

V dPt J^=i I dPt j,=i+l dPt ^^^^ 

Applying Eq. (58a) to above equation, and using the fact that the potential V{r) is inde- 
pendent of P in the limit of infinite lattice sizes and infinitesimal lattice spacing (a — > 0), 
one can obtain. 
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datV{r) \ _ 1 / datV{r) 1 datV{r) 



1^,, , /(9a a\ , , 

o^WU?-^ , (60) 



2 ' 13 

where we have used (^)^^i = || and (^)^=i = —a. Also the second term on the L.H.S. 

of Eq. (59), ( ^a^/f'' )g=i; can be evaluated as, 



V d(3t A=i 2W/5 /3 



(61) 



where the self-energy f{P)/a depends on /3 and a. Substituting Eqs. (60) and (61) into 
Eq. (59) yields 

( ^atEo{r) ^ _ ly..(da _ a\ l/dfi/3) _ fi(3)\ 

[ op, )^=,-2^^'\dp p)^2\dp P )■ ^^^^ 
Similarly, the L.H.S. of Eq. (56b) can be evaluated by applying Eq. (58b), the result is. 

From Eqs. (56a), (56b), (57a), (57b) and (62), (63) one can get the energy sum rules in 
Eqs. (30) and (31). 

Now we proceed to consider the generalization of these sum rules to the finite temperature 
case. In the limit of finite temperature (T) and infinite volume the potential V{r) would 
be a function of r and T (or /3) because of T = l/Nta{P). For example, in the confined 
phase V{r) ~ kt, with the string tension n = k.{P) for fixed lattice size Nt, as discussed in 
Sec. inc. In this case one can write, V{r) — V{r,(3). Then Eq. (60) should be rewritten 
as. 

Similarly one has 

(^)..4(^ 

The self-energy part f{f3)/a does not change in this case. Therefore, at finite temperature 
Michael sum rules in Eqs. (30) and (31) should be modified as following, 
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(67) 



These modified sum rules can be obtained from the original ones in Eqs. (30) and (31) by 
the replacement, 

We expect that the modified sum rules would account for the finite temperature effects 
on the qq system. Although Eqs. (66) and (67) are derived by considering the Wilson loop 
representation of the qq pair, we expect that they can be applied in describing the flux 
distributions /^^ in Eq. (19), which involves the Polyakov loops, because both the Wilson 
loop, W, and Polyakov loops, P{0)P^{r), represent a qq pair. 



B. Analysis of Flux Data 

We calculated the flux distributions /^^ of Eq. (19) on various lattices, which cover 
the regions of both confined and unconfined phases. The detailed analysis of the qq flux 
distributions is presented elsewhere [5] . In this part we shall compare our flux data with the 
predictions of modified Michael sum rules in Eqs. (66) and (67). 

As we mentioned in Sec. II C, we used the multihit technique [6] in the flux measurements. 
The flux data is only good when the test charge (plaquette) is away from the qq sources 
(Polyakov loops). Thus we only consider the flux data on the middle transverse slice between 
the qq pair at large separations (i.e., r > 3a). 

In the confined phase string formation is expected to occur. When the separation, r, of 
a qq pair is large enough, the flux energy stored in the flux tube per unit length between the 
qq should equal to the string tension k. In this case one can obtain the color electric and 
magnetic energy of the flux tube per unit length from Eqs. (66) and (67), 
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dEei{r) 1 fp d[aK\ , 



dr 2\a dp 

dEma{r) lfPd[aK] 



+ (70) 



of Eq. (15), e.g., ~ —8 for P ~ 2.4, which is consistent with the flux measurement 
results of Ref. [10]. The second term can be estimated from the fitting function in 



dr 2\a dp 
The corresponding flux action is 

pdjan] 
a dp 

which is the integration of the action density on the transverse plane. And the corresponding 
flux energy is 

dV(r) 

aE = del + CT^a = = n{P) , (72) 

which agrees with our expectation for string formation. 

Eq. (71) can also be obtained from the action sum rule of Eq. (47). The R.H.S. of 
Eq. (71) is, f = + P^- The flrst term can be estimated from the scaling relation 

8 In g 
ain/3 

Eq. (28), and it is negative because k, decreases with P for flxed lattice size Nt, as shown 
in Eq. (25). So Eq. (71) implies that aA ^ (Je — k for cases we consider. Also, Eqs. (69) 
and (70) show that the electric part of the flux energy cTei is positive, but the magnetic 
part (Jma negative, their magnitudes are of the same order with a^i slightly larger. These 
predictions can be summarized as below, 

(7el> and (Tina < 0, 

(7e= CTel + CTma = K > 0, 

o-A> cte] (73) 

where these flux energies (action) are deflned in the conflned phase. 

From our flux data we can calculate the values of cTei, (Jma., <^e and cx^, which are the 
energies (action) stored in the center slice per unit length between the qq pair. We flnd 
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that our data agree with the predictions of Eq. (73) quahtatively. However, it is difficult to 
measure accurately because of the strong cancellation between two terms, (7ei and (Jma, 
which have opposite signs. This problem was also discussed by Refs. [9,10]. 

In the following we proceed to study the prediction of the action sum rule in Eq. (71). 
As we mentioned above, the R.H.S. of Eq. (71) contains two terms, the first term, /tfj^, 
was predicted by the original Michael sum rules in Eqs. (30) and (31), the second one, 
is a new term, which is only predicted by the modified sum rules of Eqs. (66) and (67). We 
expect that the second term describes finite temperature effects. 

From the scaling relation a(/5) of Eq. (15) one can estimate the quantity. 



dlnp 121 ir Al\P) 

with A2^{(3) defined in Eq. (15). And by using the fitting function in Eq. (28) one can 
calculate the derivative, 

with the constants Kq and 5 given by Eqs. (14) and (29). 

By similar considerations that lead to Eq. (71), one can also obtain the prediction about 
a A from the original Michael sum rules in Eqs. (30) and (31). The result is 

I \ ^^^^ ('7R\ 

= -'^9h^' 

which is just the first term of the R.H.S. in Eq. (71). 

In Table IV we list the predictions of Eq. (71) obtained by substituting Eqs. (74) and 
(75), and the measured a a data in the confined phase. Since in the confined phase the value 
of a A should not change with the qq separation r in the limit of r — cxd, we then choose 
the a A data at a moderate value of r as the asymptotic value of a a {t — > oo), because 
error bars are large at very large r. In Table IV most a a data were calculated from the 
flux measurements of qq pair at r = 4a, in few cases we choose the data at r = 3a or 5a, 
depending on the quality of data. 
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In Fig. 7 we plot the predictions of Eqs. (71) and (76) respectively, and compare them 
with the measured a a data on lattices of Nt = 4. Prom this figure one can see that the 
measured data are consistent with the prediction of Eq. (71), but have large difference from 
the prediction of Eq. (76). Especially, in the transition region (/3 /3c) Eq. (76) predicts that 
(7a approaches zero as j3 ^ Pc- However, our data have large values in this region, which 
agree with the prediction of Eq. (71). As the temperature decreases (T — > or /3 — > 0), 
both predictions coincide. However, it is difficult to obtain data of clear signal in the small 
P region (i.e., P < 2.25). 

We also notice that in Fig. 7 the cr^ data in the confined region (/3 < Pc) have some 
discrepancies from the prediction of Eq. (71). This may be due to finite-size effects and 
finite lattice spacing effects of lattices, and to the large fluctuations in the confined phase 
because confinement corresponds to disorder. 

To compare the behaviors of cta in both phases, in Fig. 7 we also plot the cta data in 
the unconfined region (i.e., P > Pc). These data were measured at large qq separations, i.e., 
r = 6a, which still have clear signals, as shown in Table V. Since in the unconfined phase 
there is no string formation, one expects that aA vanishes at large r. From this figure one 
can see that near the transition point, i.e., P ~ 2.30, the aA data in the unconfined region 
decrease rapidly with P, and becomes very small beyond the transition region. This agrees 
with the expectation. The fact that our cr^ data in the unconfined region do not vanish may 
be due to the following factors, in the transition region finite-size effects are large, beyond 
this region the contribution from the self-energy of the qq pair still exists, because the qq 
separation, r = 6a, is not large enough to approach the asymptotic region. 

V. SUMMARY 

We studied the SU(2) finite temperature phase transition by Monte Carlo simulations. 
To transform the measurements from lattice units to physical units the lattice asymptotic 
scaling relation a{P) is extrapolated into the non-perturbative region. The behavior of string 
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tension with temperature is also studied, We find that the string tension data agree very 
well with the fitting function, k{T) = ko(1 - for T < T^. 

We also measured the flux distribution of the qq pair at various temperatures. To check 
Michael sum rules with the flux data, a complete derivation of the sum rules are presented, 
and a generalization of the sum rules is suggested to account for the flnite temperature 
effects. We found that our flux data are consistent with the prediction of the generalized 
sum rules. Our data shows explicitly that the qq flux distribution has different behaviors in 
the two phases. In the conflned phase the asymptotic value of the center slice action, a a 
(r —>■ oo), has large values. However, across the transition region, a a becomes very small in 
the unconflned phase. This agrees with our expectation that string formation occurs in the 
conflned phase, but disappears in the unconflned phase. 
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APPENDIX A: 



To prove Eqs. (58a) and (58b), we notice that Eqs. (53a) and (53b) implies that jSg, fit 
and /3, ^ are two equivalent sets of variables, that is, F{Ps:Pt) — F{Ps{P:Ot Pt{P:0)- Then 
we take the partial derivatives of F with respect to /3 and ^ respectively, one has 

(^] = (^] (^] + (^] (^] (Ai) 
(^F^ ^(^F^ ^^F^ /a^x 

Prom Eqs. (53a) and (53b) one can get that, when /3 is large (i.e., (3 — > oo). 



^ (3 + 2Nc[\^=i. (A3) 



where denotes for SU(A"), and c' the ^-derivative of c. Substituting Eq. (A3) into Eqs. (Al) 
and (A2) yields 

/dF\ /dF\ /dF\ 

(AS) 



(f) 



«=1 



||(/J + 2iVc;) + |£(-/J + 2iV4) 



F. Karsch has studied the coefficients and Cj [25]. His results show that the derivatives 
c's\^=i, cj|^=i vanish for sufficient large lattice couphng constant, /3. For example, as /3 > 2.2, 
the values of c^,|^=i and c^|^=i of SU(2) LGT are about 0.1, which is much less than the 
values of /3. For large /3 cases Eqs. (A4) and (A5) become 

fdF\ fdF\ fdF\ 

(f).rKII)..-KI£).. 

After carrying out some simple algebra, we can obtain Eqs. (58a) and (58b). 
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TABLES 

TABLE I. The correspondence of the lattice spacing a and the couphng constant /3 for SU(2) 
LGT, extracted from the string tension data of Refs. [8,19]. 



(3 




a{(3) (fm) 








2.22 




0.1981 (61) 






27.41 (84) 


2.30 




0.1616 (47) 






27.32 (79) 


2.40 




0.1210 (5) 






26.30 (11) 


2.50 




0.0843 (8) 






23.57 (22) 


TABLE IL The raw string 


tension data k measured in 


lattice units on 


lattices of the size, 


4 X 92 


X 65, 4 X ll2 X 65, 6 X 7^ 


X 65 and 6 x 11^ x 37. 






Nt = 4 


4 X 92 X 65 


4 X 11^ X 65 


Nt = 6 


6 X 72 X 65 


6 X 11^ X 37 


P 


^/Ka 




(3 


^/Ka 


^/Ka 


2.25 


0.301 (38) 


0.313 (32) 


2.30 


0.310 (48) 


0.377 (27) 


2.28 


0.242 (33) 


0.243 (29) 


2.36 


0.248 (21) 


0.284 (13) 


2.29 


0.224 (17) 


0.181 (33) 


2.40 


0.172 (30) 


0.217 (17) 


2.30 


0.216 (19) 


0.203 (17) 


2.42 


0.200 (17) 


0.177 (41) 
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TABLE III. The string tension data k measured in physical units on lattices of the size, 
4 X 9^ X 65, 4 X 11^ X 65, 6 X 7^ X 65 and 6 x 11^ x 37, which were calculated from the data 
in Table II by using the scaling relation a(/3) in Eq.( 15). 







AT A 




4X9 X DO 


4 X ii X DO 


P 






1 (i/imj 


K (LieV /imj 


K (LreV/imj 


z.zo 
2.28 
2.29 
2.30 




on 
3/. 29 

40.20 

41.22 

42.27 


1.2/ i 
1.401 
1.447 
1.494 


0.46 (i2j 
0.36 (10) 
0.33 (5) 
0.33 (6) 


O.oO (iOj 
0.37 (9) 
0.22 (8) 
0.29 (5) 






Y, = G 




G X 7*^ X G5 


G X 11^ X 37 


/? 






T (1/fm) 


K (GeV/fm) 


K (GeV/fm) 


2.30 
2.36 
2.40 
2.42 




28.18 
32.76 
36.23 
38.10 


0.996 
1.213 
1.386 
1.482 


0.68 (21) 
0.64 (11) 
0.40 (14) 
0.62 (11) 


1.00 (14) 
0.84 (8) 
0.64 (10) 
0.49 (23) 


TABLE IV. The predictions from Eq. (71) and the data of center slice action a a, which were 
measured on lattices 4 x 9^ x 65, 4 x 11^ x 65, 6 x 7^ x 65 and 6 x 11^ x 37. The quantities k(/3), /?||, 
'^Unf ^'^'^ ^s;ve the physical unit GeV/fm. The values of k(/3) were estimated from Eq. (28). 

Near the transition i")()iut (i.e.. 3,. ~ 2.30 for Ni = 1) uo sLalilo i")r(xiic(jou is obiaiuod. 








Nt = 


4 




/? 


k{(5) 


„ dlna 
'^dln/S 






(<''a)4-112.65 


2.25 
2.28 
2.29 


0.42 
0.34 
0.29 


3.03 (47) 

2.50 (48) 
2.12 (47) 


4(1) 
9(3) 
17 (6) 


7 (2) 12.18 (23) 
12 (3) 10.95 (14) 
19 (7) 10.44 (12) 


8.70 (22) 
9.97 (13) 
9.39 (10) 
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TABLE V. The a a data in the unconfined phase (/3 < /?c), which were measured on lattices, 
4 X 9^ X 65 and 4 X 11^ x 65. The data are in the physical unit Gev/fm. 



Nt = 4 

P 2.30 2.32 2.34 2.36 2.40 

(^a)4x92x65 11.77 (11) 10.40 (9) 1.78 (4) 1.28 (4) 

(^a)4xii2x65 11-03 3.81 (4) 1.77 (3) 1.32 (4) 



32 



FIGURES 

FIG. 1. Monte Carlo data for (|P|) vs. (3 with the standard Wilson action, calculated from 

lattices 4 x 9^ x 65 (circles), 4 x 11^ x 65 (squares), 6 x 7^ x 65 (triangles) and 6 x 11^ x 37 (diamonds). 

FIG. 2. Monte Carlo data of (P(O)P(r)) vs. r/a in the confined phase, calculated on the lattice, 
4 x 11^ x 65 with /S = 2.25 (squares), /? = 2.28 (triangles) and /3 = 2.29 (diamonds). These data 
are all in the confined region, with (5 < (3c {Pc ~ 2.30 for Nt = A). The plot is in a frame with 
logarithmic y axis. 

FIG. 3. Monte Carlo data of (P(O)P(r)) vs. r/a in the unconfined phase, calculated on the 
lattice 4 X 9^ X 65 with /? = 2.36 (squares), and /? = 2.40 (triangles). 

FIG. 4. The plot of k vs. T (1/fm) near the transition point Tc ~ 1.487 (1/fm). Data were 
calculated on lattices, 4 x 9^ x 65 (circles), 4 x 11^ x 65 (squares) and 6 x 7^ x 65 (triangles). The 
solid line is the fitting function in Eq. (26) with a = 0.35 zt 0.04. The string tension k is in the 
physical unit GeV/fm. 

FIG. 5. K vs. (3 for Nt = 4, the transition point is chosen as (3c = 2.2985. Data were calculated 
on lattices, 4 x 9^ x 65 (squares), 4 x 11^ x 65 (triangles). The solid line is the fitting function in 
Eq. (28) with 6 = 0.22 zt 0.03. The string tension k is in the physical unit GeV/fm. 

FIG. 6. The Wilson loop W of the temporal size na and the spatial size r. The plaquette Pi 
is outside the Wilson loop W, and the plaquette P2 is inside the Wilson loop W. 

FIG. 7. The plot of the predictions of a a vs. P from Eq. (71) (solid lines), and Eq. (76) (dashed 
lines) in the confined region {(3 < (3c)- The two solid lines represent the upper and lower limits 
predicted by Eq. (71). The data were measured on lattices, 4 x 9^ x 65 (squares) and 4 x 11^ x 65 
(triangles). For comparison the data in the unconfined region {(3 > (3c) are also shown, which were 
measured on the same lattices, 4 x 9^ x 65 (circles) and 4 x 11^ x 65 (diamonds). Here a a has the 
physical unit GeV/fm, and the transition point (3c is indicated by the up arrow. 
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